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SUMMARY 


The differential equations for the geometric layout of com- 
pression-loaded axisymmetric and cylindrical two-family filamen- 
tary structures are established. The analytical formulation is 
based upon the requirement that failure due to local instability 
occurs simultaneously in the whole structure. 

Solutions are obtained for the particular case where the body 
force due to the structure's own weight is the only load. For 
this special case, the shapes of the meridian and the cross sec- 
tion, respectively, have been determined, as well as £he pattern 
of the filaments. In addition, the weight per area covered by 
the structure is formulated and discussed. Minimum-weight 
configurations are defined. Examples for large earth— based and 
moon-based structures are presented. 
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INTRODUCTION 


Among high-perf ormance structures maximum economy of material 
is often obtained by a design in which all structural components 
fail simultaneously as the structure is subjected to a unique 
ultimate loading condition. Examples of such structures are the 
uniformly stressed Maxwell-Michell type networks, isotropically 
stressed structural " soap-film" membranes, uniformly stressed 
filamentary network structures ("isotensoids") (Refs. 1, 2), and 
the '‘Shanley tube" (Ref. 3) , which is designed to fail simultan- 
eously in a local and general instability mode. 

Analogously to the isotensoids , structures can be found 
which are uniformly compressed. These rr isocompressoids" are the 
subject of a recent study (Ref. 4). However, under compressive 
loads, failure of the structure due to instability may occur long 
before the compressive stress limit is reached. The present report 
treats this stability problem from a standpoint of buckling of the 
individual structural members. Hence, the analytical design con- 
dition imposed upon the structure is that each element buckle as 
an Euler column between network nodes upon reaching the design 
ultimate load. The term ,f isostabiloid" is used to describe the 
filamentary structural configurations satisfying this condition. 

This analysis deals with axisymmetric and cylindrical config- 
urations with axisymmetric (or uniform along the length) loadings. 
In either case no circumferential (or lengthwise) loadings are 
considered. The influences of general instability are recognized 
but only briefly discussed. 



LIST OF SYMBOLS 


A 

b 

C 

D 


d, d , 
E 



cross-sectional area of a single rod 
. 2 

" Sln P a 

buckling end-fixity coefficient 
tube diameter 

d^ dimensions of particular Sintes column (Fig. 13) 
elasticity modulus 

weight form factors for gravity- loaded isostabiloids 

stress form factors for gravity- loaded isostabiloids 

gravitational acceleration 
. 2 

= sm B 

moment of inertia of rod cross section 

curvature at apex of structure in nondimensional 
coordinates 

stiffness form factor of rod cross section 



m mass per unit planform area of gravity-loaded isostabiloids 

n number of filaments in each family 

n number of rod sections in each filament 

N number of unrestrained filament intersections 


p external meridional load per unit surface area 

*m 


p external normal load per unit surface area 

n 

P local buckling load of rod 

cr 


q 


number of meridional node lines in general failure mode 
of axisymmetric shells 


2 
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r radial coordinate 

R nondimensional radial coordinate in gravity-loaded 

ax i symmetric structure (see Eq. 5) 

s spacing of filaments at perimeter 

t w.all thickness 

x widthwise coordinate in cylindrical structures 

X nondimensional widthwise coordinate in gravity-loaded 

cylindrical structures (see Eq. 7) 

y axial coordinate in cylindrical structures 

Y nondimensional axial coordinate in gravity-loaded 

cylindrical structures (see Eq. 7) 

z vertical coordinate 

Z nondimensional vertical coordinate in gravity-loaded 

axisymmetric structures (see Eq. 5) 

Z nondimensional vertical coordinate in gravity-loaded 

cylindrical structures (see Eq. 7) 

a angle between meridian and vertical 

p angle between meridian and filament 

At length of filament between two intersections 

p density of rod material 

a stress 

Subscripts : 

. ... values at apex of cylindrical structures 

.... values at rim of structure 

o 
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BASIC RELATIONS 


Ax i symmetric Isostabiloid 

Consider an axisymmetric structure consisting of two famil- 
ies of interconnected rods as the load-carrying parts, and built 
in such a way that the arrangement pattern of one family of rod 
trajectories is the reflected image of the other family. The rod 
trajectories are referred to as "filaments", indicating that these 
may be continuous structural members, connected to other filaments 
at localized points of intersection. 

Let n be the number of filaments in each family; then the 
length of an element. At , between two intersections is, with the 
notation of Figure 1, 


At 


1 2ur 1 

2 n sing 


( 1 ) 


Assume that all rods are of the same cross section and 
material; thus the bending stiffness, El , is a constant. Let 
another constant, C , be the end-fixity coefficient for the 
single rod element. This parameter varies from one to four as 
the end fixity increases from pinned to clamped. If the rods 
are sufficiently slender, then the force in the filament causing 
fully elastic local buckling can be written as 


P 

cr 


C n El 


A-t 


or, from Eg. (1) 


C El n .2 

P = sxn R 

cr 2 w 

r 


( 2 ) 


Assume that the loading is axisymmetric. 
brium equations can be written. 


Then two equili- 


4 



Let p m and p be the external loads per unit surface 
m n 

area of the structure in the meridional and the normal direc- 
tions, respectively, as shown in Figure 1. Also, let P be the 
compressive force in the filament. Now, examine the filament 
between two intersections with filaments of the other family. 

rrr Ar 

The corresponding surface area has the size : and is 

n sma 

traversed diagonally by the filament section, A£ . Equilibrium 
requires that in the direction of the z-axis 

. , _ _ , rrr Ar rrr Ar 

A(P cosa cosB) - p : cosa + p ; sma = 0 

m n sma n n sma 

and in the radial direction 

A(P sina cosB) - P sinB — - p sina 

sma r m n sma 

ur Ar 

- p ; cosa = 0 

n n sma 

Dividing by Ar and assuming that the filaments are suf- 
ficiently closely spaced that the difference quotient can be 
replaced by the derivatives, yields 


~z (P cosa cos p ) - p — cota + p — = 0 
dr m n n 


~ (P sina cos R ) 
dr 


sinS tanB nr 

- p — 

r sma m n 


p — cota = 0 
n n 


In order to find a more convenient second equilibrium equa- 
tion, the first one is multiplied by cosa cosp , the second one 
by sina cosp , and the results are summed: 


dP 

dr 


sinfi d_ 
r dr 


(Pr sinp ) 



0 
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It is now postulated that at ultimate loading the compres- 
sive force, P , reaches the local fully elastic buckling load 
everywhere in the structure. Thus, P is replaced by P 
from Eq. (2 ) : 


cr 


C El 


2 d / 
n dr V 


sin B cosB cosa 
2 


\ TT 

) ~ P m n r 


cota + p — r = 0 
n n 


C El 


2 \ d / sin^B \ sinB d / sin~^B \ \ 
n j dr \ r 2 ) ~ r dr \ r / j 


TT COSB 

- p — r — ; — L = 0 
m n sina 


Some manipulation yields a form convenient for subsequent use 


/ 4 2 2 \ 

d I sin 3 cos B cos g \ 

dr V 4 / 


2tt p .2. . 2 

m sm P cosg cos a 


C El n' 


r sina 


2tt p . 2 

+ n sm B cosg cosa _ ^ 


d / sin^p \ _ 1_ d_ / sin u p \ 
dr \ 4 / 2 dr 2 / ' 


C EI n' 


• 6 n 
sm B 


2tt p .2 

m sm B cosB 


C EI n 


3 r sina 


[ (3) 


These two equations allow the determination of the complete 
geometry of the structure in terms of the loading and the section 
properties of the rods. 


Cylindrical Isostabiloids 

A cylindrical case, such as is shown in Figure 5, can be 
construed as the limit of an annular type of axisymmetric struc- 
ture with indefinitely increasing radius. 

The differential equations for the geometry of cylindrical 
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isostabiloids can therefore be derived directly from the equations 
for axisymmetric structures, Eqs. (3), by substituting (r^ + x) 

for r and letting r^ approach infinity. Consequently dr be- 
comes dx , and the spacing of filaments 

2tt r^ 


is introduced. The limiting process yields the basic equations 
for the cylindrical isostabiloid: 


d . . 4 Q 2 2 . 

— (sm 3 cos p cos a) - p 

dx m „ 2 


3 2 2 

s sin 3 cos3 cos a 


4 tt C El 
3 


+ P. 


n 4 tt^ C El 


sina 


2 

sin 3 cosp cosa = 0 


d 

dx 


( sm p 


cos^ 3 ) 


m 


3 

s 

2 

4tt C El 


sm 


cos3 


sma 


0 


(4) 
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GRAVITY LOADING 


Consider the structure to be subject to its own weight with 
the direction of the gravitational force in the z-direction, and 
with g as the gravitational acceleration. Let p be the density 
of the rod material and A be the constant cross-sectional area. 
Then, the loads per unit surface area in the meridional and normal 
direction are, respectively. 


P m = 
m 


p g A 


n 

rrr 


cosg 

cosp 


n 

p = - p g A 

r rrr 


sing 
cos (3 


Eqs. (3) become 


4 2 2 

d I sin p cos (3 cos g 

dr 




)- 


o P g A sin 6 cota 

Z 2 2 ~ 0 
C El n r 


d 

dr 


sin^ft 


\ _ 1 d / sin 6 g \ 

/ " r 2 dr V r 2 / " 


. 2 

9 P 2 A sm g cota _ 

Z 2 2 — U 

C El n r 


Introducing nondimensional coordinates, 

/ V 1/3 


R = 


p q a 
C El n“ 


p g 

C El n“ 


1/3 


(5) 


and defining a variable 


. 2 

h = sm g 
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yields the differential equations for the geometry of a gravity- 
loaded axisymmetric isostabiloid 


> ( 6 ) 


A parallel procedure can be applied to the cylindrical con- 
figuration. If the lengthwise coordinate is defined as y , the 
appropriate nondimensionalization is 



p q 

A 

2 

s 

„ 2 
4tt 

C 

El 

P 9 

A 

2 

s 

. 2 
4n 

C 

EI 

p q 

A 

2 

s 

. 2 
4tt 

C 

EI 


> ( 7 ) 


The resulting differential equations for the gravity- loaded cylin- 
crical isostabiloid are: 


[h^ (1-h) cos^a] - 2h cota = 0 
dX 

d [h^ (1-h)] - 2h cota = 0 
dX 


( 8 ) 
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Domes 


Equations (6) have to be solved numerically. Because closed 
dome-shaped isostabiloids are of most interest, efforts are con- 
centrated on this type of solution. The equations are clearly 
singular at R = 0 , and careful attention must be given to the 
behavior of the solution at this point. The differential equation 
of the meridian is 


az 

dR 


cota 


If domes exist, that is, if a = 90° for R = 0 , it must be 
possible to develop the meridian, Z , and consequently cota, as 
well as h in power series in R: 


cota = ^ ^ ^ R 1 
i=0 


CO 



j=0 


By applying the differential equations (6) to these series 
it can be shown that there is a unique solution 


and that 


cota = X R + 0(R ) 


h = h 2 R 2 + 0 (R 4 ) 



The meridional shape near the center, therefore, becomes 


Z = 


Z 

o 
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Thus represents the nondimens ional curvature of the meridian 

at the apex. Call it k . 

For the purpose of numerical integration the differential 
equations (6) are rewritten as 

2 . 3 

da _ h cota - R 

dR “ Rh(l-h) 


dh 

dR 


h(2-h) + R^ cota 
R( 1 - | h) 


and the first term approximations 


cota = k R 

. 1 2 
h =-R 

Z = Z + | R 2 


are used to start the process. 

In the course of integration the independent variable is 
switched from R to h in order to avoid a division by zero 
when h reaches the value of 2/3. This value is significant 

because the meridian ends there 555 °)- Thus, on a 

closed-dome axisymmetric isostabiloid, the filament angle, p , 
can vary only between 0° and 54.74°. 

Because of the closed-dome specification, the filament path, 
p(R) # and the meridional shape, a(R) , consequently also Z(R) 
and the height- to-semispan ratio, Z Q / r 0 / are completely defined 

if the nondimens ional curvature at the apex, k , and the termina- 
tion of the integration, R , a , or p , are selected. 

o o o 

For example, Figure 2 shows side and top views of a dome- 
shaped isostabiloid, specified by the curvature, k = 2 , and 
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the filament angle at the rim, = 54.74°. 

The nondimens ional curvature at the apex, k , is an inde- 
pendent design parameter. It is needed for the numerical evalua- 
tion, but it provides little information about the structural 
configuration. For the selection of a particular dome, therefore, 
a more practical design parameter, the height- to-semispan ratio, 

z /r , is plotted in Figure 3 versus the integration limit, R 
o o o 

with k as parameter. In addition, the values of a Q and 8 0 

corresponding to each combination of R q and k are indicated 

by lines of constant value. Therefore Figure 3 is a design chart 

for dome-shaped isostabiloid configurations because now any 

desired combination of the parameters, a ,8 , R , z /r , and 

o -o o o o 

k , can be picked to select a unique configuration. 

As an example, domes with a height- to-semispan ratio, z ^/ t q 

= 1, 2, and 3, and 8 Q = 54.74° are seen to be characterized by 

k-values of approximately 3, 6, and 10, respectively. Their 
meridians have been computed and are shown in a normalized form 
in Figure 4. 


Tunnels 

For the cylindrical configuration, attention is directed to 
singly-connected tunnel-type shapes. Equations (8) can be re- 
written as 


iL_ ^ h 2 (1-h) sin 2 aj = 0 
dX 




2h cota = 0 


The first differential equation can be integrated to yield 


. 2 
sm a 


b 2 (1-b) 
h 2 (1-h) 


12 



2 

where b = sin P at the apex, where a = tt/2 . Substitution 

a 

into the second equation yields 


or 


u( h 2 ( i-h>) + ;*- h ; 


d 

dX 


dX = + 


-1 = 0 


b (1-b) 


tyUb (i- |h) 


Vh 2 Cl- 


dh 


(1-h) - b (1-b) 


If the integration starts at h = b , tunnels closed at the 
apex are produced because the jt sign in the differential equation 
implies a symmetry with respect to the origin of the X-coordinate. 

The cross-sectional shape of the cylinder and the lengthwise 
pattern of the filaments are determined by the two equations 


and 


dZ = cota dX = 


± (1 - f h) dh 


dY = 


tanl ^ 

sina 


= + 


h 3/2 (1 - | h) 




dh 


h) - b (1-b) 


The polynomial under_the radical in the denominator of the form- 
ulas for dX and dY can be written as 


h 2 ( 1-h) - b 2 ( 1-b) 


(a-h) (h-b) (h-c) 


where 
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and 


c 


1-b 

2 



) 


Now the coordinates, X and Y , can be written in terms of ellip- 
tic integrals 


X 


Y 


b 




(h-b) 


(h-c) 



dh 


while Z = J (1 - ~ h) dh 

b 

Note that a and b are exchangeable, and that further- __ 
more a = b for b = 2/3 . But this, too, is the point where X 
and Z reverse their directions in the integration. Therefore, 
again, at h = 2/3 the rim is reached. For values of b 
smaller than 2/3 , the filament angle, (3 , increases from its 

value at the apex, |3 , to 54.74°, while for b larger than 2/3, 

0 decreases from the apex to the rim. 

An investigation of the nondimens ional curvature at the 
apex, k , shows that 


k 


d 2 Z 

dX 2 


d (cota) 
dh 


dh 

dX 


1 

b (1-b) 


and therefore can never be smaller than 4. 
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As in the dome structure, the geometry of the tunnel is 

determined by specifying the behavior at the apex (it is P , 

b , or k whichever form is preferred) and the end of the 

intergration (a , P , or X ) . For example. Figure 5 shows 
o o o 

an isometric projection of a cylindrical isostabiloid specified 
by the two values of the filament angle, P: 30° at the apex 

(equivalent to b = 0.25 or k = 5.333) and 54.74° at the rim. 

Figures 6 and 7 present the design charts for the configura- 
tion of cylindrical isostabiloids, with filament angles at the 
apex, (3 , smaller than 54.74° in Figure 6 and larger than 

54.74° in Figure 7. Again, from this chart the design parameters 
are chosen which pertain to cross-sections of tunnels with a 
height-to-semispan ratio, z /x , of 1, 2, and 3, shown in 
Figure 8. 
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MASS OF THE GRAVITY -LOAD ED STRUCTURES 


Because gravity is the only external load, the mass of the 

structure can be determined from the forces exerted upon the edge 

support. Each filament contributes to the total vertical load 

the amount, P cosa cos 8 , where P is defined by Eq. (2) . 

cr o o cr ^ 

Thus the average structural mass per unit planform area, m , is, 

for domes, 


m 


2n P cosa cosB 
cr o o 


2 

g n r o 


and for tunnels, 


m = 


2 P cosa cosg 
cr o o 

s g x 


Let the moment of inertia of the individual structural 
member be related to its area by the formula 


I = K_ A 
f 


where is a shape factor. Also replace P^ r by Eq. (2), then 

for domes , 


2n 3 C E K A 2 2 

m = sin (3 cos|3 cosa 

4 o o o 

g it r 


(9) 


and for tunnels, 


8rr 2 C E K f A 2 2 

m = sin p cosP cosa 

3 o o o 

g s x 


( 10 ) 
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I 


Now, from Eg. (5), for domes, 

SL 


A = 


C E K f n 2 \ R o 


(ID 


Also, from Eg. (7), for tunnels. 


A = 


2 x 

p g s o 


2 X 

4tt C E K f o 


(12) 


If these expressions for A are substituted into the appro- 
priate formulas for mass, and the spacing at the rim of the dome 


2tt r 

o 

s = 

o n 


is introduced, the following results are obtained. 


For domes, 


and for tunnels, 


p g s r 
K o o „ 

m = — f 

2 D 

tt C E K„ 


2 

p g s x o 
m = — 

tt C E K f 


(13) 


(14) 


The functions, f and f , are expressed in terms related to 

the structural configuration alone and, therfore, are called 
weight form factors: 
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2 

sin 8 cosf3 cosa 
_o __o o 

D " „ 6 


2 

sin (3 cosP cosa 
o o __o 

f T — 2 

2X 


They are plotted in Figure 9 versus the height- to-semispan 

ratio z /r or z /x , with the nondimensional radius of the 
o o o o __ 

dome, R , and the nondimensional semispan of the tunnel, X 
respectively, as parameter (full lines) . For easier interpreta- 
tion, lines of constant k-values are also indicated (dashed) in 
the two graphs. 

It has been found above that for constant k-values the 
filament angle at the rim, P Q , increases or decreases steadily 
with growing R Q — and X - values. Consequently, Figure 9 shows 
that for constant values of k the weight form factors of gravity- 
loaded axisymmetric isostabiloids become a minimum when P 0 
reaches its maximum possible value of 54.74°. The global minimum 
of the weight form factor, f =6.07 , is achieved for an iso- 
stabiloid dome with a height-to-semispan ratio z 0 /r Q of approxi- 
mately 0.53. A view of this "minimum weight" isostabiloid is 
shown in Figure 10. 

In the case of cylindrical structures, only values of f 
for 8 smaller than 54.74° are shown because for larger values 
the factor, f T , is considerably larger. Note that unlike f , 
the minimum of f for constant k-values is not always coupled 
with the largest possible filament angle at the rim. The global 
minimum, f = 3.6 , however, does belong to a configuration 

min 

with an angle, £ , of 54.74° and a z /x of 0.6. 

o o o 

Examination of Eqs. (13) and (14) shows that the operative 
parameter for material efficiency is E/p 2 . 

Of particular interest is the fact that the structural mass 
is proportional to the spacing of the filaments. This implies 
a weightwise advantage of finely textured structures. However, 
a limit is provided by the lowest order failure mode of general 
instability. 
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DESIGN LIMITATIONS 


Recall the intrinsic assumptions of the foregoing analytical 
design procedure, i.e., 

- The instability failure takes place in the fully linear, 
elastic range, and 

- The instability failure mode of the structure which 
characterizes the design consists of nodes at each filament 
intersection. 

Hence, for the procedure to yield a valid design it is necessary 
that 


- The stress at ultimate loading must remain at or below the 
proportional limit of the structural material, and at or 
below the local crippling stress of the structural section 
employed, and 

- No "general" instability failure modes occur at loading 
levels below that which corresponds to the postulated design 
failure mode. 

The two design limitations resulting from these requirements will 
be discussed separately below. 


Stress Limits 

The stress limitation of isostabiloid design is subject to 
rigorous analysis. It can be shown that the maximum stress is 
reached at the rim. Using Eqs. (2) and (11), and (2) and (12) 
yields for domes: 

ct = p g r f (15) 

°D ° D 

and for tunnels 

°o = P g X o f T (16) 

T 
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The functions f D and f T are stress form factors because 

they depend, again, on the_conf iguration alone. They are plotted 
in Figure 11. Note that f is presented only for values of 

0 smaller than 54.74° because, as in the case of the weight form 
a _ 

factor, f for 3^ larger than 54.74° assumes considerably 
higher values. 

No distinct global minimum for f^ and f appears in 

Figure 11. Indeed, a closer examination shows that for decreasing 

parameter values of R or X , as_well as for increasing k- 

values the corresponding minima of f approach a value of 1.5 

at z /r or z /x equal to 0.75. But, as it may be checked 
o o o o 

in the design charts. Figures 3 and 6, in the same process the 
filament angle at the rim, 3 q , becomes very small and the 

corresponding isostabiloids theoretically degenerate into domes 
with plane, non-intersecting radial members and into tunnels 
consisting of plane arches perpendicular to the tunnel axis. 

In addition. Figure 9 shows tremendously increasing weight form 
factors for the corresponding configurations. Consequently, in 
cases where the stress limit really might be reached, a compro- 
mise between weight and stress consideration has to be sought. 

A look at the stress equations provides an idea when the 
stress becomes pertinent: the semispan of isostabiloid structures 

will be stress limited at 
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1 


(17) 


r 

o 

max 


a 

max 



and 


x 


max 


a 

max 1 



(18) 


a 

Note that is the specific working stress, and that 1/f 

P 9" 

assumes a value of 2/3 or less. Hence, the largest, stress 
limited structure will have a span of approximately 4/3 times 
the specific working stress. 

For example, an earth based aluminum construction 
(pgearth = 0-1 lb/in. ^) operating at a maximum safe stress 

level of 10,000 psi would yield a span of 130,000 in. or 

approximately two miles 1 On the moon surface (pg = 0.0165 

2 moon 

lb/in. ) a similar structure would be stress limited at a span 

of almost 13 miles. Even with only half the values of l/f , 

assuring weight form factors close to the global minimum, the 

stress limited structures assume enormous proportions. 


General Instability Limits 

For the cylindrical structure, it is evident that the lowest 
instability mode will tend to have nodes only at the foundations 
and at the apex. Thus, the only " isostabilo id " would be the 
extremely degenerate case where the filaments do not intersect 
at locations other than at the apex; auxiliary braces will be 
required to prevent premature general instability failure. Since 
the dead weight loading associated with these braces has not been 
considered in the analysis, the design procedure will at best 
yield a first approximation to the desired optimum structure. 

For domes, the situation is more complex. Let n be the 
number of rod sections contained in each filament between apex 
and rim. Then for a structure with 2n filaments, the number 
of unrestrained filament intersections is 
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N = n (n - 1) + 1 


The elastic properties of the structure thus can be expressed by 
a stiffness matrix of order 3N or 6N , depending on the nature 
of rotational constraints at the intersections. Eigenvector 
extraction yields the buckling modes for the corresponding 3N 
(or 6N ) critical loads. The associated computational problem 
for any but the simplest structures is enormous and a rigorous 
parametric treatment is beyond the scope of this study. However, 
one simpler way to determine approximately the conditions under 
which general failure occurs is explored in the appendix. 
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APPLICATIONS 


Two sample designs are given to demonstrate the applications 
of the isostabiloid design procedure. 

Consider earth based isostabiloid domes. One dome shall have 
the optimal proportions for minimum structural weight, shown in 
Figure 10, the other shall be a design characterized by a low 
stress form factor and the same height-to-semispan ratio as the 
minimum weight dome. Its configuration is shown in Figure 12. 

Both domes shall be constructed from aluminum alloy tubing. An 
ultimate design factor of 4 shall be applied to the gravity 
loading, and the design shall be balanced to produce local crip- 
pling of the tubing at the materials proportional limit, , 

in the most severely stressed rim members. At the same limit 
stress, Euler buckling of all members is to occur. The size 
and detail construction of domes meeting these conditions shall 
be determined. 

The pertinent non-dimensional parameters are tabulated below: 


TABLE I 


Minimum Weight 


R 

o 

Z 

o 

3 

o 


a 

o 



Conf igurat ion 

0.59 

0.31 

54.7° 

47.1° 

6.07 

3.21 


w-Stress-Factor 

Configuration 

0.39 

0.20 

18.3° 

42.3° 

20.2 

1.68 
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The materials properties are assumed to be 

E = 10 7 psi 

pg = 0.101 lb/in. 3 
E 


a — 40 , 000 psi 
P 

To meet the required design factor of 4, the maximum per- 
missible working stress under simple gravity load is 


a = 10, 000 psi 
o 

The size of the dome is now determined by Eq. (17) 


r 

o 


o 

o 


P9 t 


f 


D 


This yields 

r^ = 2,570 ft for the minimum weight configuration, and 

r = 4,890 ft for the low-stress-factor configuration, 
o 

The diameter, D , and wall thickness, t , of the tubing is 
chosen such that local crippling of the wall will occur at the 
proportional limit. Hence, from Reference 3: 


a = a = 0.4 E t/D = 40,000 psi 
cr p 


or 


t/D =0.01 


and 



3.98 
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The structure shall be built such that its members will 
buckle at a load equal to 4 times that which is produced by its 
own gravity. Equation (5) can be rewritten to yield the diameter 
of the structural tubing which will satisfy this condition; 



. 3 

r o 4 ( pg) 

3 2 

R CEn 
o 


and, assuming C — Is 


D = ^ ^ for the minimum weight configuration, and 
n 

D = 1 for the low-stress-factor configuration, 

n 

The weight of the structure is found from Eq. (13), 
which can be rewritten, by substituting 

2tt r 


g = 4 g E , and s q 


n 


o 


which yields: 


2 8 C “V 2f D 

w “ mg E r = " — ^!k; — 


or 

W = X - -^ lb for the minimum weight configuration, and 

n 


47 9 x 10^ 

W = 2 lb for the low-stress-factor configuration. 

n 
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Table II and III below summarize the properties of isosta- 
biloid domes for various values of n , 


TABLE II 

MINIMUM WEIGHT CONFIGURATION 


n 


36 

72 


100 

250 

500 

1000 

D 

(in.) 

187.5 

93 

.75 

67.5 

27 

13.5 

6.75 

t 

(in. ) 

1.87 5 

0 . 

9375 

0.67 5 

0.27 

0.135 

0.0675 

W 

(10 6 lb) 

313.5 

157 

113 

45.2 

22.6 

11.3 





TABLE : 

III 






LOW-STRESS 

-FACTOR 

CONFIGURATION 




n 

24 

36 

60 

100 

180 

360 720 



D (in.) 

1400 

937. 

5 560 

336 

187 

93.75 46.9 



t (in.) 

14 

9.375 5.6 

3.36 

1.87 

0.9375 0.47 



W (10 9 lb) 

20 

13.3 

7.98 

4.79 

2.66 

1.33 0.665 


Moon Structure: 

Future bases located on the surface of the moon will require 
shelters for personnel, material and equipment. An isostabiloid 
dome may be employed as a scaffold in the initial phase of the 
construction of such shelters. Since the structure will be 
transported from the surface of the earth, minimum weight and 
storage volume will be required. 
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As an example of such a structure, consider a "minimum weight" 
isostabiloid dome with a diameter of 333 feet (r Q = 2000") and 
apex height of 88 feet. It is to be designed with a load factor 
of 3 for operation in the moons gravitational field of 1/6 g E , 
and the dome is to be made from flattenable and coilable 
"Sintes-Tubes " (Ref. 5) . 

The cross section of this type of structural member is shown 
in Figure 13. As material of construction 17-7 PH Stainless 
Steel is selected with the following properties: 

6 3 

E = 28 x 10 psi, (pg ) = 0.28 lb/in. , a = 180,000 psi 

E y 

The maximum working stress at the rim is given by Eq. (15) 

Q o = r o 6 (Pg E } ? D = 300 Psi 


This structure will obviously not be stress limited. 


The wall thickness/convolute diameter ratio of the structural 
tubes, t/d , is determined by the condition for flattening: 
Allowing a maximum strain of .33% during the flattening process 
requires d/t = 300 and yields a stiffness form factor = 15.6. 

The (earth) weight of the structure is derived from Eq. (13) 


E 


" mg E = 


4 2 

r o (PV 

n C E K_ 


4360 

n 


lb 


and the cross sectional dimension of the tubing (see Fig. 13) 
become 


. 49.8 . 

d, = in. 

1 n 


, 45.3 . 

d = in. 

2 n 


t = 


0.0798 

n 


in. 
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Reasonable structural sections are obtained by choosing 
n = 72 . The structural properties and dimensions for this 
design are given in Table IV below: 


TABLE IV 


333 FT DIAMETER MOON DOME STRUCTURE, n = 72 


Design Factor 

3 

Weight (lbm) 

60.5 

d x (in.) 

0.69 

d 2 (in.) 

0.63 

t (mils) 

1.1 


28 


CONCLUDING REMARKS 


A mathematical formulation suitable for analytical design 
of structural networks with uniform local instability margins 
has yielded solutions for specific, idealized cases of simple, 
two-family "axisymmetric" and "cylindrical" networks. Such 
structural configurations may form the basis for the design lay- 
out of practical structures. Not considered were the effects of 
unavoidable local and general imperfections upon the optimal 
choice of design configurations, and the effect of multiple 
loading conditions. In addition, only one type of loading was 
treated in detail, and the influence of general instability 
phenomena was only shortly discussed. A fertile field for further 
theoretical and experimental studies is indicated. 



P + dP 



Figure 1. Isostabiloids - Geometry and Forces of an Element 
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APPENDIX 


GENERAL INSTABILITY OF DOME-SHAPED ISOSTABILOIDS 


The lowest general instability mode of a dome-shaped iso- 
stabiloid may be determined approximately by replacing the network 
structure with an "equivalent" uniform shell of revolution. Solu- 
tions to shell buckling problems of this type are available in the 
published literature . 

For a dome-shaped shell of uniform wall thickness, t s , and 
subject to axial loading, the lowest instability failure mode 
seems to be inextensional , according to Reference 6, and it will 
generally exhibit meridional node lines (see inset of Figure A-l) . 
The number of node lines, q , depends, among other parameters, on 
the wall thickness-to-semispan ratio of the shell, t s /r Q . Fig- 
ure A-l, derived from data of Reference 6, shows this dependence 
for the particular case of a paraboloidal dome with a meridional 
angle at the shell rim, a Q , of 45°. 

Now, by defining an "equivalent" wall thickness for isosta- 
biloids the data of Reference 6 may be applied to these gridwork 
structures provided the meridional shapes of the shell and the 
isostabiloid are similar. 

From inspection of the buckling mode it appears that the 
bending stiffness in circumferential direction is the dominant 
elastic shell property controlling the buckling mechanism. While 
the bending stiffness per unit transverse width of the isotropic 
shell is uniform 


Et 


El ' = 
s 


12 ( 1-v ) 


(A-l) 


it varies with location and direction for an interconnected grid- 
work of rods, forming an isostabiloid, and becomes in the circum- 
ferential direction 


n . 3 

El 1 = El sm f 

c ur 


(tanp+J^otp) 


(A-2) 
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El and GJ are the bending stiffness and the torsional stiffness 
of a single straight rod. 


The equivalent wall thickness of an isostabiloid, t e g is 
now defined by combining equations (A-l) and (A-2) and integrating 
over the complete structure to obtain an average value 

r 

o 

I ' dr 
c 

o 


t 3 = 12(l-v 2 ) — 
eq r. 


/ 


Using Eq. (5) to substitute r Q the equivalent wall thickness is 
found in its final general form as 




. 4 GJ . 2 
sin p+— sin p 

R cosp 


dR 


(A-3) 


For the special case where the structural members are made 
from circular tubing Eq. (A-3) reduces to 

R q 

12 / pgl 2 a \ 1//3 1 f ( 1 — y 2 ) sin 4 g+(l-v) sin 2 p , 

TT y 1 CE / R Q / R cosg 

o 

which, using v = 0.3 , yields for the minimum weight configuration 




mw 


0.231 


TT ^ CE J 


Referring to Figure A-l, it can now be postulated that if 
q ^ 2n then the corresponding isostabiloid will indeed fail in 
the assumed manner. If q ^ n it may be suspected that the low- 
est instability mode of the corresponding isostabiloid involves 
motions of joints normal to the surface and that the isostabiloid 
design will not be valid, unless additional constraints (such as 
braces, stay wires, etc.) are used to prevent such failure modes. 
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The foregoing approximations must be used with caution since 
they do not account for unsymmetrical instability modes that may 
be caused by unavoidable geometrical imperfections and by side 
load components; nor do they account for failure modes that are 
associated with tangential strains. 
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